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Abstract 

We report some new exact instantons in general relativity. These 
solutions are Kahler and fall into the symmetry classes of Bianchi types 
VIo and VIIq, with matter content of a stiff fluid. The qualitative 
behaviour of the solutions is presented, and we compare it to the 
known results of the corresponding self-dual Bianchi solutions. We 
also give axisymmetric Bianchi VIIq solutions with an electromagnetic 
field. 



1 Introduction 

Solutions to field equations in imaginary time have been of interest for many 
years and early work formalised the relevance of such systems to tunnelling 
problems 0. Much of the current interest in instantons dates from the 
work of Belavin, Polyakov, Shvarts, and Tyupkin (BPST) [Q], where an exact 
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solution to the SU(2) Yang Mills equations in imaginary time was found. The 
relevance of these solutions was then clarified by t'Hooft [Q and by Jackiw 
and Rebbi ||^ who showed that the BPST instanton solution is a mediator 
for the transition between different vacua of the gauge theory. 

The BPST instanton was found by looking for solutions with (anti) self- 
dual field strength. This led to the suggestion that analogous solutions may 
exist in gravity 0, where the duality governs the curvature two-form. The 
usefulness of such a duality is that it not only makes the curvature automat- 
ically solve Einstein's vacuum equations 0, but also reduces the equations 
to first order. A geometry which possesses a symmetry also simplifies the 
problem. An important class of symmetries is the Bianchi classification of 
spatially homogeneous 3-geometries|Q. Some solutions of this type with the 
property of self duality are known and can be found in reference ; we shall 
extend the collection of known solutions by finding metrics with the symme- 
try of Bianchi types VIq and VIIq that satisfy Einstein's equations with a 
stiff fluid, that is where the pressure equals the energy density. This means 
that we cannot use the simplification of self duality, however we do use it as 
a guide to find an initial ansatz. We shall also see that this ansatz allows a 
natural complex structure such that the metric is Kahler. 

The importance of gravitational instantons emerges in the path-integral 
formalism of quantum gravity. In this approach one considers a weighted 
sum over all Euclidean metrics, where the weight is the exponential of the 
Euclidean action. Instantons, as solutions to the classical field equations, 
are expected to provide the dominant contribution to such a sum, and this 
allows a saddle-point approximation to be used in the neighbourhood of these 
solutions. 

The organisation of the paper is as follows. First we define the system of 
equations to be solved, and describe the self-dual vacuum solution. We solve 
the Bianchi VIIq Einstein equations in the presence of a stiff fluid. Next, we 
follow Misner |T^ by writing the action in the familiar form of a particle in 
a potential, in order to describe the qualitative nature of our solution. The 
Bianchi VIq solution follows from the Bianchi VIIq results. 
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2 The Bianchi VIIq equations 



The Bianchi universes are four-manifolds with embedded three geometries 
which are spatially homogeneous but anisotropic; they each possess a three- 
dimensional, simply-transitive, isometry group To exploit this isometry 
we write the metric using the left-invariant one- forms of the symmetry group, 
a"; as we are interested in instantons, we use a Riemannian signature, 

= df2 + a\f){a^f + h\f){a^f + c\f){a^f. (1) 

The left-invariant one-forms satisfy the structure equations, 

d^'^ = \C%,^' A (2) 

For the case of Bianchi VIIq, the structure constants take the values 

C32 = = C\.^ = —C\i = 1. (3) 

These metrics support Killing vectors (right-invariant vector fields) with the 
same structure constants as the Euclidean group in two dimensions. The 
self-dual solution will be of use to us in constructing an ansatz, so we shall 
calculate the torsion-free connection forms, B"^ 0, assuming they are metric 
connections, Qab = — ©ba-Thus, we rewrite the metric in terms of a new time 
variable df = abcdr and introduce an orthonormal basis, 

ds^ = a^{T)b^{T)c\T)dT^ + a\T){a'f + b\T){a'y + c''{T){a^)^ (4) 
= a{T)b{T)c{T)dr, uj^ = a(r)a\ 00^ = b(r)(T^ 00^ = c{T)a\ 

We find the following non-vanishing connection forms: 

^0 = -a'/{abc)uj\ e\ = -{a^ + b^)/{2abc)uj^ 

e\ = -(3'/{abc)uo\ e\ = (a' - b'') / {2abc)uo^ (5) 

^0 = -i/{abc)uo\ e\ = -{a^ -b^)/{2abc)uj\ 

where we have defined a = ln(a), (3 = ln(6), 7 = ln(c) and ' refers to differ- 
entiation with respect to r. We now assume the connection to be self dual, 
which in turn implies that the curvature is self dual so 601 = 623 plus 
permutations. 

a' = -l{a'-b% P' = l{a'-b% y = ^(a^ + &^). (6) 
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These equations show that a(r)6(r) is a constant, and give a metric, 



ds' = — sinh(2r) [dr' + {a-'Y) + coth(r)(ai)^ + tanh(r)(a^)^ (7) 

where A is an integration constant. 

We now derive a non-vacuum solution, (which requires that self duality 
is lost). However, we use the self-dual solution to lead us to the ansatz of 
a(r) = l/6(r). In order to show that this solution does retain the Kahler 
property, the metric is rewritten by introducing the complex forms C, defined 
by 

= u^ + iu^ (8) 
= oj^ + ioj'^, (9) 

where the tu" are those defined in (^). The metric now takes the form 

ds2 = CV + C^C^. (10) 

The Kahler form is then found to be 

= IQalC^C' (11) 

= Aabc^dr A + 4abcr^ A (12) 

The closure of the Kahler form, d^K = 0, then requires (ab)' = 0, so our 
ansatz leads to a Kahler manifold. 

We now want to consider the full Bianchi VIIq equations, without requir- 
ing self duality. These can be found in refs. [Til], [1^ where we are now using 
a Lorentzian signature, 

ds' = a\t)b\ty{t)dt^ (13) 

(\n{a^))" + - b" = (14) 
(inib^))" - + b" = (15) 
(ln(c2))--(a2-62)2 = (16) 

(iniabcf)" = 4((lnc)-(lna6)- + (lna)-(ln6)-) + 2e(a6c)l(17) 
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where e > is the energy density, and • denotes differentiation with respect to 
t. Now we take the hint from the self-dual solution and assume a{t) = l/b{t), 
but we do not take e to be zero. Then equations (|T^,|T3|,[TB|) become, 

= -^^Hy) (18) 

= sinh^(y/2), (19) 

where we have used x = 2t, = exp(?//2) and = exp{h). Equation (p!^) 
can be solved to give, 

cosh{y{x)) = 1 + 2—cn'^{x/k' + a \ A;^), (20) 

where k is the modulus of the elliptic function, k' is the complementary 
modulus, and a is an arbitrary constant. Our notation conforms to that of 
With (|0l) substituted into (^) we find 

= ^cn^(xA' + a|^^). (21) 

We can integrate this equation once by noting that the elliptic integral E 
satisfies, [|14| , 

Ij^E Um{x/k' + a I k^)) = 1 + cn'ix/k' + a \ A;^), (22) 

allowing us to write (^I]) as 

^ = ^e( am(x/k' + a I P)) - x + d. (23) 
ax k' ^ ' 

A change of variables to ^(x), and the introduction of a new function h{x) 
defined by, 

h{x) = hi{x) + J dx{-x + Ci), ^ = am{x/k' + a \ k'^) (24) 

means that we have 
dhi 



E{i\k')l^l-kHm\0 (25) 
h,{i) = F'{ak')E{k)/{2K{k)) + ln{em\k')]/e[0]). (26) 
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The integration has been done using equation (630.02) of |T^. This fur- 
ther simphfies when we reahse F(^ | k"^) = x/k' + a. The conservation of the 
energy-momentum tensor, T'^''.^^ = 0, gives the energy density e as 



e = {(3/{ahc)f = {(3/c)\ (27) 

where /3 is a constant. The fourth Einstein equation (p!7|) reduces to 

h = - (iyf/A + 4/?2) . (28) 

But we aheady know y{t) is from (pOD, and the left-hand side comes from 
(PTD, so we have the two consistency relations 

+ = 1, e/k'^ = (29) 

In what follows, we shall use the value < /5 < 1, with k = 1/3/ y/1 — (3'^ and 
k' = 1/ \/l — (3'^. Now that we have a relation between j3 and k we can rewrite 
(pOD using some elliptic function identities. We also make use of the fact that 
a is an arbitrary integration constant which just shifts the origin of the 
solution. It proves convenient to choose a = y/1 — P'^K{P), turning (|20|) into 
(see appendix 0) 

cosh(y(t)) = 1 - 2/5W(2t I (3^). (30) 

This is where we make a connection with the instanton solution. It is seen 
that if we take t to be an imaginary parameter, t = ir, then the right-hand 
side of (30) becomes greater than or equal to one, meaning that ?/(r) is real. 

cosh (?/(r)) = 1 + 2(3V{2t \ 1 - (3^)/cn^{2T \ 1 - (3^). (31) 

To proceed then, we are required to do an analytic continuation on (^) from 
t to r (see appendix 0). 

h{r) = 2t\1 - E{(3)/Km - vrrV {k{(3)K{^1 - Z^^)) (32) 
+ ln[e(2r I 1 - P^)] + Cir + 

This completes the solution for the Bianchi VIIq instanton with a stiff fluid. 
The solution is defined in the range K{^J1 — P"^) < 2r < K{^/Y^^]3^), with 
c(r) going to zero at these limits. At these points then we see that the volume 
of the homogeneous 3-surfaces goes to zero, defining the end-points of the 
instanton. It can also be seen that the curvature of the instanton diverges at 
these points. 
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3 Qualitative nature of the solution 



Let us now examine the system of equations (p!4| - |T7|) in order to obtain an 
effective action. Using imaginary time, r, we see that these equations can be 
derived by extremising the action, 

Ses = / dr (a'/3' + /?V + «V ' ' b'f) • (33) 

Again, we are using a = ln(a), (3 = ln(6), 7 = ln(c) and ' for differentiation 
with respect to r. The form of this action coincides with the Einstein-Hilbert 
action for the Bianchi VIIq metric, found by calculating the curvature of 
the connection forms (|^) so that this effective action is proportional to the 
gravitational action, as expected, 

-^grav = / v/^-Rd^x + boundary terms. (34) 

The boundary terms are designed to cancel the second derivatives coming 
from the Ricci scalar term, |jl5|, ||16|, in order that the equations of motion 
derive from the variational principle even on a boundary. The next step is 
to put the effective action (|33|) into a familiar form, so we can visualise out 
solution. This is achieved by using the following Misner variables |T7[ , 

a(r) = exp(fi(r) +/3+(r) + V3/3_(r)) 

6(r) = exp(fi(r)+/5+(r)- V3/3_(r)) (35) 

c(r) = exp(^](r)-2/?+(r)) 

We then find that the effective action, in terms of the new functions, takes 
the familiar form for a particle moving in a potential well, (^) 

-3 J dr + p'J + ^ exp(4fi) exp(4/5+) sinh2(2V3/3_)^ (36) 

The P± functions evolve to create the time-dependent potential 
V(fi,/5+,/9_) = exp(4n)W(/5+,/5_). Fig. |I] illustrates the (^-independent part 
of the potential, with a ridge along the /3+ axis. The behaviour of the two 
solutions outlined above, the self-dual and stiff fluid instantons, are given 
in Fig. 1^. The figure shows the self-dual solution (paths Ii and I2) evolves 
along the ridge and then falls off at some point determined by the arbitrary 
parameter A of the solution, (|^). The solution for the stiff fluid contains three 
parameters, k, Ci and C2.The general behaviour is shown in Fig. which 
shows the solution evolving toward and then rolling over the ridge. 



'cff 
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4 Bianchi VIq solution 



The Bianchi VIq structure constants differ from Bianchi VIIq (|^) by a single 
minus sign: 

C32 = = — C'ls = C'si = 1- (37) 

This changes the symmetry group from the two dimensional Euclidean group, 
E(2), to E(l,l), and also changes one of the equations of motion, ([T6|). Using 
the same metric as (P, but with the left-invariant one forms of E(l,l), we 
find the new equation, 

(ln(c2)) - - (a2-62)2 = 0. (38) 
Again we make the ansatz a{t) = l/b{t) giving 



dx^ 



1 + sinh^(?//2), (39) 



where we have again taken x = 2t, a? = exp{y/2) and = exp{h). Compar- 
ing equations (|l^) and (|39D , we see that the solution for h{x) in the Bianchi 
VIq case is the same as for the Bianchi VIIq case, except that there is an 
addition of . However, the Bianchi VIq solution differs in other respects. 
By looking at the fourth Einstein equation, (p^), we find a different relation 



between k, k' and the parameter (3 of (p7|). Using (pHj) and the new h(t) we 
may derive, 

k'^ + e = i, kyk'^ = -{i + (3^). (40) 

Using methods similar to the previous analysis one may then derive. 



[3 
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cosh(|/(r)) = l + 2(l + /32)cs2(^2^1+/32r|^^j (41) 
Mr) = -2(l-i5;(/3/v/r+^)/ir(/?/v/r+^))r2 (42) 



8 



This instanton is defined in the range < r < E{l3 / j ^/T^^W ■ Fol- 

lowing the procedure outhned for the Bianchi VIIq case, we may also derive 
the self-dual solution, 

ds^ = ^ sin(2r) (dr^ + {a^f) + coi{T){a^f + tan(r)(a2)2^ (43) 

The Hamiltonian method used to describe the qualitative behaviour may also 
be adapted. On doing so we find that the solution is described by the effective 
action as (p6|), only now the sinh^(2A/3/?_) gets replaced by cosh^(2-\/3/3_), 
using the same allocations as (|35D. The effective potential for the /3± is thus 
qualitatively the same as Fig. |l|. The paths that this solution describes in 
the (3± plane are shown in Fig. ^, where the examples of self-dual solutions 
are labelled Ji and /2, with a, h and c being stiff fluid examples. 



5 Bianchi VIIq solution with electromagnetic 
field 

In this section we give an analytic solution of Bianchi type VIIq in the pres- 
ence of an electromagnetic two-form satisfying Maxwell's relations. Using 
the orthonormal basis for the metric we write Maxwell's relations as, 

dF = 0, d*F = (44) 

which are valid here because there are no charged sources. We consider here 
a two-form field strength of the form 

F = Ecu° A cu^ + Buj^ A cj^ (45) 

representing an electric and a magnetic field in the '3' direction. One quickly 
sees then that (0) leads to both E and B being constant for the Bianchi VIIq 
metric. When it comes to the field equations for the metric coefficients we 
have only a minor change. The right-hand side of ( p^ - |T6D no longer vanishes 
but becomes r]Uc^, where 77 is -|-1 for (0), ([T5|) and —1 for (p^ . The constant 
U that appears is given hyU = {EY + {BY (see |T^). There is a clear choice 
of ansatz for these equations, namely the biaxial ansatz with a = b. Taking 
this ansatz leads to 

(Ina^)" = _t/c2 (46) 
(Inc')" = Uc\ (47) 
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and with the definition = exp(^^) the equation for c gives, 



e" = Uexp{e), (48) 

taking note that we are now using Euchdean time, r = it. We can trivially 
perform the first integral of (|48| ) by multiplying both sides by 6' giving. 



d^ 



/D^ + 2Uexp{e) 
Integrating this and using (|^) we find 



dr. (49) 



c\r) = ^^^^ (50) 

2U 

a\T) = exp(Cir + C2)^^— ^ (51) 

l-exp(±D(r + ro))' ^ ^ 

with Ci, C2, D and tq all constants of integration. The 4th Einstein field 



equation (eq.(15) in |T8|) gives the constraint Ci = D. The qualitative be- 
haviour of this electromagnetic Bianchi VIIq model follows from the previous 
discussion. The extra terms rjUc^ in the equations of motion correspond to 
the addition of an extra term \Uc^ in the effective action (133|), which has the 



effect of tilting the potential in ( pG] ) in the (3+ direction. The trajectory in 
(/3+, (3-) space is clear in this biaxial example, /5_ = so the solution rolls 
along the axis. 



6 Conclusions 

Exact solutions of the Einstein equations are only possible in the presence 
of symmetries or structural specialisations. When the metric signature is 
Riemannian it is possible to find vacuum solutions when self duality is com- 
bined with three-dimensional homogeneity. We have used the structure of 
these solutions to find a number of non-vacuum solutions with stiff matter 
and electromagnetic field sources when self duality is necessarily relaxed. 
These solutions are of Bianchi types VI and VII. 
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A Appendix 

In order to transform fEOf) to f|31| 



cosh(2/) = 1 - 2/?^cn^ 2t/k' + a 



(53) 



\ 1 - V 

we use the transformation of elliptic functions for negative parameter, as 



found in [13 



cosh(y) = 1 - 2p^cd^ ((2t/k' + a)/y^l-/?2 I . (54) 

We use the fact that a is arbitrary and define, — [3'^ = K{P), and find 

cosh(|/) = 1 - 2pV{2t I p^). (55) 

The final step is to perform the continuation t = ir which gives us (|HT|). 
Here we continue equation (pSj) to imaginary time t = ir, 



h{t) = {2t/k' + afE{k)/{2K{k)) + In {Q[2t/k' + a \ A;^]/e[0]) - 2^ + dt + ^2.(56) 
We use the definition given by Cayley of the theta function, [p^ , 
0(m I 



exp 



0(0 I F) 
^ h{t) 



Define the double integral in (|57|) to be, 

r2t/k'+a 

h ■ 



1 /"^ 

-n2(l - E{k)/K{k)) dz dysn\y \ k^ 

2 Jo Jo 



-{2t/k' + af -21" + Cit + C2 

l-2t/k'+a rz 

—k"^ / dz dy sn^{y \ k"^ 

Jo Jo 



(57) 
(58) 



dz dy sn (y \ k ). 
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(59) 
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Change t to ir and make two changes of variables, z = 2iv /k' + a and y = /k' + a 
to obtain 



A(r)^Ar + g-4A"/-dv/;d^ ^^,^^^^\_^,^ . (60) 

In deriving this we have used the relation between fc, k' and (3 (p9D; note 
that k is imaginary. The double integral can be rewritten in terms of theta 
functions by making the following change of variables, ^ = iu/2 and v = iy/2. 
This leaves us with a double integral which may be related to a theta function 
with imaginary argument. This can be expressed in terms of a theta function 
with real argument using a result in ||19[, and hence we obtain (0). 
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Figure 1: U{(3+,(5J) 
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^ 

Figure 2: Various Bianchi VIIq instanton solutions in the /3+ — /9_ plane. 
The Ii_ 2 are examples of self dual solutions, with curves a, b, c being three 
stiff fluid instantons. 
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